Influence of Low Frequency Noise on Macroscopic Quantum Tunneling in 

Superconducting Circuits 
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The influence of low to moderate frequency environments on Macroscopic Quantum Tunneling 
(MQT) in superconducting circuits is studied within the ImF approach to evaluate tunneling rates. 
Particular attention is paid to two model environments, namely, a pure sluggish bath and a slug- 
gish bath with additional 1/f noise. General findings are applied to Zener flip tunneling, a MQT 
phenomenon recently predicted and observed in a superconducting circuit implementing a quantum 
bit. 
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I. INTRODUCTION 

The decay of a metastable system through macro- 
scopic quantum tunneling (MQT) has been studied in 
view of conceptual questions in quantum theory already 
in the 1980s Recently, it has regained new inter- 

est e.g. in the context of quantum information processing 
in solid state based systems 0, Q or as an effective de- 
tection mechanism in shot noise measurements @. In 
corresponding electrical circuits, current-biased Joseph- 
son junctions are used as building blocks where the phase 
difference across a junction is the only relevant degree of 
freedom. The switching from the zero voltage state can 
then be visualized as the dynamics of a fictitious particle 
moving in a cubic potential such that at low temperatures 
the escape is dominated by quantum tunneling through 
the potential barrier. The crucial impact of the electro- 
magnetic environment on this process has already been 
elucidated twenty years ago and then led to the devel- 
opment of the " standard" theory of quantum dissipative 
systems [g. 

To date, the most powerful approach to calculate tun- 
neling rates in dissipative systems is the so-called lmF- 
method which relates the escape rate to the imaginary 
part of the free energy of an unstable system 0,0- Basi- 
cally, it can be seen as the generalization of the method 
to extract the lifetime of an unstable state from an imagi- 
nary part of its resonance energy Q . The theoretical pre- 
dictions of the ImF theory including dissipation have 
been thoroughly tested in the 1980s in Josephson junc- 
tion circuits Q|. In these studies the dominant effect of 
friction was due to dynamical environmental modes with 
ohmic spectral density. These modes generate at very 
low temperatures T a rate enhancement proportional to 
T 2 as compared to the zero temperature limit [l(J, UJ$ ■ 

In recent experiments with superconducting circuits 
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for the implementation of quantum bits, the role of 
low frequency modes in the environmental spectrum has 
turned out to be substantial E! El Eil d Bath 
modes that are inert on the typical experimental time 
scales or show only slow to moderate dynamics basi- 
cally determine the dephasing time for oscillations of co- 
herent superpositions of qubit states. In many cases, 
the corresponding noise spectrum displays a 1/f depen- 
dence over a broad frequency range attributed e.g. to 
low frequency bistable charge fluctuators or electromag- 
netic fluctuations in control lines. In fact, 1// noise has 
been known for years to be always present in mesoscopic 
devices 0, 0. While detailed theoretical studies on 
its influence on the decoherence process in qubit devices 
have been given, apart from some qualitative estimates, 
much less is known about its impact on MQT. 

In this paper, we analyse MQT processes subject to 
noise from very low up to moderate frequencies. One may 
expect that quantum fluctuations in the bath are thus 
suppressed. For two model environments, namely, a pure 
sluggish bath and an environment with 1/f characteris- 
tic, we show in detail to what extent this is actually true 
and how such a situation is incorporated into the ImF 
theory. These findings are further applied to a particu- 
lar kind of MQT phenomenon which has been predicted 
E3 and observed [2(J recently in the so-called quantro- 
nium circuit Q|. This system, a superconducting circuit 
implementing a two-level system (qubit) with a readout 
by MQT, realizes an extension of the standard MQT sce- 
nario: It describes a fictitious particle with a spin-i as an 
internal degree of freedom such that the Zeeman splitting 
of its levels is position dependent; in certain ranges of pa- 
rameter space the Zeeman levels cross under the barrier 
leading to spin flips while tunneling. The MQT rate is 
then substantially enhanced since effectively the particle 
has to penetrate a smaller barrier. 

The article is organized as follows. In Sec.|n]we give a 
brief account of how tunneling rates in dissipative quan- 
tum systems are calculated within the ImF approach. 
Then, for the generic case of a cubic potential the two 
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models of low frequency environments are analysed, a 
sluggish bath model in Sec. II I II and a model that incor- 
porates additional dynamical modes with a 1// charac- 
teristic in Sec. II VI The MQT readout in the quantronium 
device is addressed in Sec. [V] particularly the Zener flip 
MQT. At the end some conclusions are given. 



calculation of the free energy F — (—1/(3) \ogZ and thus 
of the partition function Z of an unstable system. Most 
conveniently, this is done within the path integral repre- 
sentation which gives, after tracing out the bath degrees 
of freedom exactly, the partition function of the reduced 
systems as 



II. FREE ENERGY METHOD FOR 
TUNNELING RATES 

The stochastic motion of a one-dimensional degree of 
freedom q with mass M in a potential field V(q) is de- 
termined by a classical Langevin equation 

f* dV 
Mq + M ds 7 (t-s)q(s) + —=Z(t). (1) 



The damping kernel j(t) and the Gaussian fluctuating 
force with = are related by the fluctuation- 

dissipation theorem (£(i)£(s)} = Mk B Tj(t - s). The 
situation where a particle is initially confined in a 
metastable well and, as time elapses, leaves it due to ther- 
mal activation has been studied in a variety of systems 
0, l2l| . In the low temperature range below a certain 
crossover temperature To, however, quantum tunneling 
is known to be the dominant escape process. Particu- 
larly, for an undamped system one has Tq = hu>o/{2irkB) 
where u>q denotes the frequency for oscillations around 
the minimum of the well (plasma frequency). 

To describe quantum barrier penetration in presence 
of a dissipative environment, one starts from a sys- 
tem+reservoir model H = Hq + Hb + Hj, where 



H 
H B 
Hi 



^ 2m 



C .; .1 ; 



(2) 



This model of a bath consisting of harmonic oscillators 
captures linear dissipation and equivalently mimics an 
environment with Gaussian statistics. In the limit of a 
quasi-continuum of oscillators the effective influence of 
the heat bath is determined by the temperature T and 
the spectral bath density J(w). Classically, for this model 
one regains the Langevin equation where the Laplace 
transform of the damping kernel 



dtj(t)e~ 







is related to the spectral density via 

1 f°° duj J{uj) 2z 
M J ~ 



70) 



LO LO 2 + Z 2 



(3) 



(4) 



To obtain the tunneling rate, the usual procedure is 
to apply the so-called Im-F method. It is based on the 



Vq e 



-S[q]/n 



(5) 



q(-hf3/2)=q(hf3/2) 



with the effective Euclidean action functional S[q] 
So[q] + Si[q] where 



So[q] = dr 
J-hp/2 



f <? 2 + v® 



(6) 



represents the bare system and 



Sl[q] 



np/2 ph/3/2 

dr I dr' K(T-T')q(T)q(T') (7) 



-hp/2 



-hp/2 



originates from the coupling to the heat bath (influence 
functional). The influence kernel follows from 



Kir) 



dw cosh[uj(hp/2-T)} 
- J{U) sinhW/2) (8) 



and is directly related to the analytic continuation to 
imaginary times of the force-force autocorrelation func- 
tion of the bath Lit) = (£(*X(0)) where £ = J2i c i x i' 1 - e - 
K(t) = L(-ir)/h. Further, due to the relation Q the 
influence kernel can be inferred from the classical damp- 
ing kernel ja, [2^| . 

The above path integral 10 sums over all closed paths 
in the imaginary time interval [— Uj3/2, h/3/2]. In case 
of sufficiently high potential barriers a semiclassical ap- 
proximation applies so that the partition function is dom- 
inated by all periodic minimal action paths and Gaussian 
fluctuations around them. It turns out that while fluctu- 
ations around a well minimum of a metastable potential 
are stable, those around a barrier top are not. As shown 
by Langer 0, due to an analytic continuation this gives 
rise to an imaginary contribution to the partition func- 
tion Z w Z m + iZ\y = Z w (l + iZ-b/Z w ) which, even though 
it is exponentially small compared to the well contribu- 
tion Z w , must be taken into account to finally obtain the 



escape 



rate [H Q 



r 



-InxF 



2_Z± 

h(3 z w 



(9) 



For the non-dissipative case and T = the method 
is completely equivalent to the WKB formalism and in 
other ranges its results have been verified by full dynam- 
ical approaches (see also @). Its theoretical predic- 
tions have been thoroughly tested in comparison with ex- 
perimental data over broad temperatures ranges and for 
various systems . The advantage of the method is that 
it enables us to calculate the decay rate of a metastable 
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system from a purely thermodynamic quantity without 
considering the complicated real time dynamics. 
The archetypical form of a metastable potential is 



V(q) = yw ? 



1 - — 



(10) 



with a well located around q = and a barrier top at 
q h = 2<7o/3. The barrier height is V& = 2/27Mwq9o an( i 
the well frequency cJo- Metastability justifying a semi- 
classical approximation to the partition function then 
means that V& 3> hu>o, k^T . For a current-biased Joseph- 
son junction the known tilted washboard potential for its 
phase takes locally the above form of a cubic surface. 

Accordingly, in the low temperature domain Z\> is dom- 
inated by the so-called bounce orbit, a minimal action 
path that runs in the time interval [—h/3/2, h/3/2] through 
the inverted barrier potential —V(q) with minimal action 
Sb. The corresponding tunneling rate turns out to be 



r = 




(ii) 



where Dq and D[ are functional determinants captur- 
ing Gaussian fluctuations around the constant well orbit 
q(r) = and the bounce, respectively. In the latter one, 
the contribution from a zero mode direction in function 
space is omitted. In absence of a heat bath and for T = 
the above expression leads to 



r = 6,„y / |e- 36 "/ 5 



where 



v = 3_ 

is the dimensionless barrier height. 



(12) 



(13) 



III. SLUGGISH BATH 

While the impact of moderate to fast noise on tun- 
neling rates has been elucidated analytically and nu- 
merically in the past |(|, here, we focus on the low 
frequency limit and start with a sluggish bath. This 
means that the fastest modes available in the environ- 
ment are still slow on the thermal time scale and com- 
pared to the system's well frequency. Specifically, we 
assume a spectral density J < (u>) with a cut-off frequency 
u) c , i.e. J<(oS) — 0, for u> > u> c , obeying hf3uj c <C 1 and 
Lu c /uJo <C 1. As a consequence, the influence kernel be- 
comes approximately constant 



Tlj3 J TT 



dto J<(u>) _ 2 



and the influence functional simply reads 



Si[q] 



h/3 



h/3/2 



h/3/2 



dr q(r) 



(14) 



(15) 



Any bath coupling establishes correlations with corre- 
sponding memory times that relate the time evolution of 
the system at the present with its past. The time inde- 
pendent influence kernel of (|15fl can be understood as the 
limiting case where the frequencies of the environmental 
degrees of freedom are are so low, that the critical corre- 
lation time ranges over the whole interval [— h/3/2, h/3/2]. 

The sluggish bath can now most conveniently be 
treated by exploiting a Hubbard-Stratonovich transfor- 
mation and introducing an auxiliary variable. Accord- 
ingly, the partition function is obtained as 



Z(a„ 



dx exp 



2ira x J_ oa \ 2a,, 



Z(x), (16) 



where Z(x) is the path integral 



Z(x) 



Vq e 



-sM/n 



The action functional S x [q] = So[q] + S x j[q] with 

rh/3/2 
-h/3/2 



rh/3/2 

S x ,i[q] = X I drq(r) 
J -h/3/2 



(17) 



(18) 



can be seen as a linear coupling between the " centroid" of 
the imaginary time orbit and a constant "external force" . 
The total partition function is thus represented as an av- 
erage over individual partition functions where the con- 
tribution of each one is weighted according to a Gaussian 
distribution with width a x — 2J S / f3. 

Now, by a proper coordinate shift for a fixed value of 
the auxiliary variable q = q + Q(x), the S x j contribution 
in the effective action is absorbed in a potential contri- 
bution of an effective bare system. While the measure of 
the path integral remains unchanged, the parameters of 
the shifted potential become n-dependent, namely, 



V(q,x) 



M 



wo(aOV 1 



q~o{%) 



(19) 



where q and ujq are defined by demanding that 
xq + V(q) = V(q, x) + V(Q(x)) + xQ(x). With the ab- 
breviation 



/ Qx 
a ( x ) = tM + TTFX 



(20) 



one obtains the coordinate shift and the modified param- 
eters as 



Q(x) = |[l-a(x)] 



ujq(x) = uj oy /a(x) 
qo(x) = q a(x). 



(21) 

(22) 
(23) 



The partition function for fixed x now takes the form 
Z(x) = e -W«W)+*M <fvqe-i^ x \ (24) 
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where So[q; x] denotes the action of the bare system © 
with the x dependent potential V from (|19fl . 

Now, in line with the semiclassical approximation, for 
sufficiently small width a x (details see below) the average 
over the auxiliary variable can be treated together with 
the fluctuations around the bounce q = qs + Sq in steep- 
est decent approximation : The full action is expanded 
around the bounce qb for x — up to terms of second 
order. This results in 

S [q;x] = Sb(x — 0) + 4— x + 3 „ x 2 

UJq MUIq 



1 

+ 2 



d rSq^[-d 2 T + ^ B '- X = 0) 



M 
~2 

-ft/3/2 v 
+ 0(x 3 ,Sq 3 ). 



M 




where Sb is the minimal action corresponding to qb- In 
this second order approximation the last term in (|26|l is 
independent of x and the remaining path integral over 
the fluctuations yields the known functional determinant 
for x = 0. The total rate then factorizes: The x = 0-part 
coincides with the non-dissipative tunneling rate T(/3) in 
a cubic potential; the remaining x-dependent terms in 
(|26[1 enter into the average over the Gaussian distribu- 
tion. Accordingly, the total rate reads 



T s (j,6)=T(6)I s (j,6), 



(26) 



where we introduced 
9 = 2ir/(uJohl3) and j = 
rection factor 



dimcnsionlcss quantities 
2J s j (Mwq), and the cor- 



(i-i + !r) 



(27) 



The correction factor provides a restriction on a x 
which in turn defines the range where the steepest de- 
scent approximation for the x-average is justified. From 
j — Qj8/2TT <C 1 one finds 



Js 



Muil 



< 1 



(28) 



For example, in case of low frequency Drude friction 
J<{u>) = M^lulo 2 /(lo 2 + lu 2 ) this leads to ^lu c /lo 2 <C 2 
which due to uj c /ujq -C 1 allows even for stronger friction. 
Upon further inspection, one sees that sufficiently below 
the crossover temperature 8^1 the full rate l|2t)|l dis- 
plays a dominating linear temperature increase due to its 
correction factor, while the bare rate approaches its zero 
temperature limit Tq (|12|l since deviations become expo- 
nentially small in 2tt/9. This linear enhancement in T is 
a specific feature of a sluggish bath u) c Tif3 <C 1. We note 
that thus in (|27l) the limit T — ► cannot be performed 
for fixed w c , but must include also ui c — > 0. Accordingly, 
zero-point fluctuations of the bath modes, which lead to 
the well-known counter term and renormalize the bare 
potential , do not appear here. 



o 

3 




0/0, 



o 



FIG. 1: MQT rates in presence of a sluggish bath vs. tem- 
perature for various values of the coupling constant j — 
0.0,0.1,0.3,0.5 (from top to bottom). Shown are results ac- 
cording to 1201 1 (solid), numerical results (dashed) and per- 
turbative results according to 1291 (dotted). 



An alternative procedure to calculate the tunneling 
in the limit of very weak dissipation and temperatures 
T To relies on a perturbative treatment around the 
non-dissipative bounce for T = 0. The first order cor- 
rection to the bare bounce action Sb is determined by 
the influence functional Si [q] in (|15fl evaluated along the 
T = 0-bouncc qb- Here, we obtain 



Si[q b ]/h 



408 



2vr 



(29) 



in accordance with the result obtained above in l|26() up 
to terms linear in j. Hence, the expression l|27|l applies 
also to larger values for j beyond the validity of the first 
order perturbation theory (see below). In this context we 
mention that the same result can be obtained by consid- 
ering fluctuations in a control parameter of the rate, e.g. 
the external bias current for Josephson junctions, and 
then taking the average over the change of the dominant 
exponential factor with an appropriate Gaussian distri- 
bution pl| . 

In order to illustrate the accuracy of the expression 
(|26|l . we evaluated the rate numerically according to the 
general formula in the case of weaker dissipation and 
a sluggish bath (see fig. 0). The bounce is gained from 
its Fourier coefficients q b (r) = ^2 Qn exp(if„r) according 
to the equation of motion in the inverted potential (see 
[9| for details). For the comparison, we further approx- 
imated T(6) k r(0) in l|26l) which is correct up to ex- 
ponentially small corrections for low temperatures. The 
agreement with the analytical result is excellent even for 
j = 0.5 and sufficiently low temperatures, while the per- 
turbative result displays negligible deviations from the 
numerical data only for j < 0.1. Obviously, the dissipa- 
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lb 



FIG. 2: Probability for escape vs. bias current for a Joseph- 
son junction with dimensionless barrier height v = 3 and 
various values of the coupling to a sluggish bath j — 
0.0 (solid), 0.1 (dashed), 0.3 (dotted). 



tive rate is always enhanced roughly linearly compared 
to the bare one due to the impact of the sluggish bath 
fluctuations. 

In experiments on Josephson junction circuits typically 
the escape rate is not measured directly, but rather the 
probability for escape as the height ib of an external bias 
current pulse varies, while keeping the width At in time 
of the pulse constant, i.e., 

-r(i 6 )At 



P{k) = 1 



(30) 



The parameters of the cubic potential l|10|) depend for 
ib close to the critical current io of the junction on the 
bias current via = u>(i b = 0)[1 — (ib/io) 2 ] 1 ^ 4 an d 

V b (i b ) = V b (i b = 0)(1 - i b /io) 3/2 - In fig. H the corre- 
sponding " s-curves" are depicted for various strengths of 
the coupling. Since in presence of slow bath modes the 
rate is larger than the bare one, the s-curves are shifted 
towards smaller values of the bias current and exhibit a 
smaller slope. 



IV. SLUGGISH BATH AND 1//-D AMPING 



annoying) source of noise ^^|. While in principle one 
could start with a truly microscopic model as in |T3, |24| , 
here, we proceed with a somewhat simplified procedure 
and effectively mimic the 1// noise by a proper spectral 
density of a heat bath. We write 



J(w) = J<(w) + Ji//(w) 



(31) 



where J< with J< (cj) =0 for u) > oj c comprises all slug- 
gish environmental modes and Ji// with Ji/f(uj) = 
for u) < uj- and ui > w+ contains its dynamical modes. 
For convenience, we introduced here in addition to oj c 
the lower cut-off o>_ with u c < uj- to have better con- 
trol over the two domains of the spectrum. Clearly, the 
separation into a sluggish part and a dynamical one is 
somewhat arbitrary and depends on the specific experi- 
mental situation, i.e. on the relevant time scales involved 
and the measurement protocol imposed. We will see at 
the end, however, that our main findings are independent 
of these latter cut-off frequencies as long as they lie in a 
range oj c , w_ <C ujq. The upper cut-off oj+ is taken to 
be only somewhat smaller than u>o such that 2ir/uj+h{3 
is still sufficiently larger than 1. For instance, for typi- 
cal experimental values ujo ~ 50 GHz and T ~ 20 mK 
corresponding to uqTi(3 « 20, one may assume w+ to lie 
in the GHz range as well, i.e. ui + ~ 1 GHz, leading to 
2Ti/Lu + h(3 w 15. 

Now, the power spectrum of bath fluctuations follows 
from the symmetrized force-force auto-correlation func- 
tion of £ = c a x a as 



1 

2^ 



(32) 



Experimentally, background charge fluctuations in meso- 
scopic circuits have shown to give rise to a behavior 
S^(u>) = sq/uj over a broad frequency range (1// noise). 
According to (|31() the spectral density J\if of the dynam- 
ical modes can then be calculated using the fluctuation- 
dissipation-theorem 



so 



9(w — w_)6(w + — lu) — hcoth(u}hj3)Ji/f(uj). (33) 



As a consequence of (|3*T1) the influence kernel K{t) 
splits into two contributions where the usual procedure 
is to split off from the second one a part local in time, 



Having analyzed the influence of the coupling to slug- 
gish bath modes we now turn to the case where in ad- 
dition to a pronounced low frequency part a tail of slow 
to moderate frequency modes are present. In particular, 
we consider dynamical fluctuations the power spectrum 
of which display a 1// characteristic. 



A. Spectral density and influence kernel 

In present experiments on superconducting qubits 1// 
noise has been found to be the most dominant (and most 



where 



and 



2 7 

K(T) = j^-k(T) + J r .S( T ): (34) 



Jf = I 2 J i//( w ) 



7T LU 



oo 

E 

n=—oo 



v n \ 7(K|) e i; 



(35) 



(36) 
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Here, : 5(t) : is the periodically continued 5 function 
with period Uj3 and v n — 2im/%j3 denote the Matsubara 
frequencies. 

The Fourier coefficients of the dynamical part of the 
influence kernel arc obtained from the spectral density 
via the relation ©. This way, for a spectral density 
with dynamical modes one finds from l|33[l ^1 //(<*>) = 
Ki/f(uj) M 0(o> — lu-)0(uj + — lo) where in the relevant 
frequency domain Ki/f(bj) — sot&nh(ujh(3) / (hu>M) w k 
is basically constant in frequency and taken to be inde- 
pendent of temperature as well. Consequently, the am- 
plitude of the power spectrum sq increases linearly with 
temperature T in agreement with an electron trapping 
mechanism for 1// noise 0|- Eventually, one arrives at 



7(W) = 



dog 



1 + 



(37) 



In particular, |^ ra |7(|^n|) = for n = 0, while for n =/= 
due to v n /oj- 3> v n /u+ 3> 1 one expands z-VilTd^nl) ~ 
(2k/ n) log(w + /cj_) — (k/-k)(lu^ —u) % _)/v\. The constant 
factor (2k /it) log(w + /w_) leads in the time domain to a 
contribution proportional to [: 5(t) : — l/h/3] where the 
part containing the 5 function exactly cancels in 134fl the 
contribution proportional to J/ . In turn, this means that 
a spectral density corresponding to 1// noise does not 
generate a contribution in the influence functional local 
in time and renormalizing the potential. 

Now, the influence functional Q takes the form 



Sj[q] 




h/3/2 



ft/3/2 



drq(r) 



(38) 



h/3/2 i-H/3/2 

dr / dr' q(r)k(T - T r )q(T') 
-H0/2 J -h/3/2 



with 



Js = J a + (Mk/tt) \0g(dJ + /uJ-) 



(39) 



being an effective sluggish bath coupling and a non-local 
damping kernel 



k(r) = C 



M 

hp 



n— — oo 



l 



where in the sum the n = term is omitted and 

C = k(w+ - JL)/ir « kuj\/it. (41) 

In fact, the above Fourier series can be summed up ex- 
actly and yields k(r) — J2t 



h(r) 



MC, 
2hf3 
M(hf3 
24 



ki(r) with 

h/3 l2 



(21 + 1) 



for t 6 (lh/3, (I + l)Ti0) 



(42) 



This way, the influence of the sluggish modes is described 
by one effective parameter J s independent of the specific 
cut-off frequencies u> c and w_, and the only relevant fre- 
quency scale of the bath is lo + . 



B. Tunneling rate 

With the influence kernel at hand, the tunneling rate 
can now be determined. The sluggish modes leading 
to the static contribution in l|39|l are treated analogous 
to the procedure described in the previous Sect. II I II 
Note that the corresponding constant coordinate shift 
does not change the form of the dynamical part since 
/ drk(r) = 0. The effective action is obtained as 
S [q;x] + §i[q\ with 

*fr0/2 rfrp/2 

dr / dr' k(r - r') q(T)q(r') (43) 

'-ft/3/2 J-%0/2 

with the x-dependent potential given in (jl9|) . 

For the case of weak dissipation and low tempera- 
tures damping due to dynamical modes is treated per- 
turbatively to first order in the dimensionless coupling 
k = k/wq by substituting the undamped T = 0-bounce 
into the influence functional <Sj. Clearly, due to the in- 
dependence of the bare action, this orbit becomes also 
a;-dependent and reads 

qo(x) 



S[q; x] 
Si[q] 



Qb,x(r) 



cosh 2 (a)o(a;)T/2) 



(44) 



Then, for the total partition function one arrives at 



Z(<r x ,a,0) 



1 



2n<j r 



dx e 2 <*x 



X e -S'' S 'i//,dyn(x) Z ^ 



(45) 



where Z(x) is specified in (|24l) and S , 1 /j dyn (x) = Si[qb. x }- 
Again the x-integration can be carried out by expanding 
the exponent up to second order in x. The resulting rate 
expression is given by 

r s ,i//(j,M) =r(0) i sA ,f(l~K,e) (46) 

with the dimensionless temperature 9 — 2ir/ujohf3 and 
the dimensionless couplings j — 2J s /Mu)q and k — k/uJq. 
The correction factor is now obtained as 



(40) I.,i/ f (J,K,d) 



1-3 



x cxp 



/ 1 - 3 + SjO/2'K 
276 jv 2 + 



cxp [-S 1/fAyn (0)/h] 

//,dy„(°)/90 



1 - j + QjO/2-K 



(47) 



Here dyn denotes the derivative of the damping term 
with respect to x and is of first order in k as well as 
Si/f, dyn- Apparently, the influence of the environment 
leads to two competing effects: The static modes give 
rise to a rate enhancement, while the contribution of the 
dynamical ones depresses the rate. To get further insight, 
we thus evaluate the action Si//,dyn(0) analytically by 
exploiting the results in l|43(l and (|44|l . and obtain 



S'l//,d y n(0)/S 



Mkqlu\ ( ai 



2a, 



2uj hf3" 



(48) 



7 



with constants a\ — 26.319..., 0,2 = 4.00.... In the low 
temperature range ivqTi(3 3> 1 this reduces to 

Si/f, dy n(0)/n « v— ^±u; + hp (49) 

7T Wo 

where we used (|41fl . An analogue calculation for the 
derivative of the damping term yields 
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^l//,dyn(0) 

3 



G 



Muj^q 



S'l// i dyii(0) 

bi 



(50) 



52.6379. 



2b 2 



&9 = 6.00 and 



with the constants 61 = 
63 = 2.00. 

We can now compare the leading order terms in the 
two exponents of the correction factor (|47|l . With re- 
spect to the static part, it is sufficient to consider only 
the dominant term 548jv/ir and there to focus on the 
contribution stemming from Ji/f, namely, Si/f^t = 
1 089vn log( lu + /u- ) / (71- 2 w 2 ) . Accordingly, 



£l//,stat _ 481og(w + /w_) 



Si 



//,dyn 



»i, 



(51) 



which reveals that the dynamical contribution of the en- 
vironment is basically negligible compared to the static 
one. Hence, the system-bath coupling parameters obey 
the relation j jk ^> 2(w+?i/3/27r) 2 . For the tunneling pro- 
cess a low to moderate frequency bath with a 1// char- 
acteristic can thus be treated as static and always leads 
to a rate enhancement which roughly grows linearly with 
temperature. The physical reason for this is simply that 
the environment has only a time period of the order of 
1/loq to probe the system while tunneling. On this time 
scale the dominant part of the bath modes is basically 
frozen, while the contribution of the dynamical ones re- 
mains small for a 1// spectrum. A changeover from a 
short to a long time behavior cannot be observed in con- 
trast to the decay of coherent superpositions in super- 
conducting qubits [TIITIITIIT^I. 

In order to illustrate these findings we show in fig. [3] 
the ratio T s l /f(6)/To for varying temperature and dif- 
ferent damping strength j = k. The analytical expres- 
sion obtained in (|4^|) coincides within this range with 
the numerical data. As seen in the inset, the influence 
of the dynamical modes to suppress the MQT rate is 
very small and only increases slightly at very low tem- 
peratures. Note again, that for fixed lo c the limit T — > 
cannot be reached. Typical temperatures in experiments, 
however, lie in the range about 20 mK corresponding to 
0/0o ^ 0.1... 0.3. 



V. QUANTRONIUM QUBIT 

The " quantronium" is a solid state based qubit setup 
consisting of a Cooper pair box whose Josephson junc- 
tion is split in two small junctions with Josephson energy 



&d 6h- 

~a 4 



0.3 0.4 




FIG. 3: Rate enhancement due to a sluggish bath with 1/f 
noise vs. temperature. The couplings constants are j = k = 
0.05 (solid), 0.1 (dotted), 0.2 (dashed), and w+/w = 1/50. 
The inset shows the suppression only due to dynamical modes 
(j = 0), but with the same values for k 



Ej, delimiting an island with capacitance C and charging 
energy Ec = (2e) 2 /2C The two lowest lying states 
of the box with energies Eg and E\ encode a qubit and 
in the circuit give rise to loop currents of opposite direc- 
tion. In parallel with the box is a third big Josephson 
junction (with E'j ^> Ej and E' c <C Ec) which serves as 
a detector. For the read-out this big junction is biased by 
an external current pulse so that the total bias current 
seen by the junction depends on the state of the qubit. 
The measurement consists of adiabatically driving the big 
junction into the regime where MQT sets in. Due to the 
exponential sensitivity of the tunneling rate on the total 
bias current, the two qubits states can efficiently be dis- 
criminated. In the quantronium two dominant sources of 
noise influence the qubit, namely, charge noise and phase 
noise. The latter one is mainly due to noise in the big 
read-out junction and there, originates e.g. from fluctu- 
ations of the external bias current. In the sequel, we are 
interested in this low frequency noise affecting the phase 
of the big junction during the tunneling process. 

The analysis of the device becomes particularly trans- 
parent in the charging regime (Ec 3> Ej) where the two 
qubit states are determined by superpositions of zero or 
one excess Cooper pairs on the island. Then, the Hamil- 
tonian of the box-subsystem can be written in terms of 
Pauli matrices and by measuring all energies in units 
of E'j the total dimensionless Hamiltonian of the circuit 
reads 



h = ec a z —ej cos 



a x + ^-cos(d)-i b 6 (52) 



with the 
(E c /E'j)(N g 



dimensionless parameters 
- 1/2), ej = Ej/E'j, and i b ■■ 



Here, N g is the reduced gate charge and I b 



ec = 
= HI b /2eE'j. 
denotes the 



external bias current. 4> represents the reduced magnetic 
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flux in units of h/2e and M — C'E'j/ie 2 is the mass 
of the artificial particle of the read-out junction with 
momentum P and conjugate phase 9 the dynamics 
of which takes place in a tilted washboard potential 
— cos(0) — ibO. For further details we refer to |J, Il9j . 
Due to the coupling between the two subsystems, the 
qubit and the big junction, the dynamics of 9 can be 
seen as that of a particle with two internal states. 

Now, we consider the situation when the qubit is ini- 
tially (if, = 0) prepared in its ground state. By rising the 
bias current if, > adiabatically the phase 9 can be seen 
as a classical variable with negligible kinetic energy due 
to the large capacitance C . Depending on the spin state, 
the artificial particle then evolves on adiabatic potential 
surfaces A± (9) obtained by simply diagonalizing h in spin 
space for M — > oo |19(. When if, is close to 1, however, 
the particle may tunnel out of the potential well. In this 
case the Hamiltonian is most conveniently represented in 
the spin basis at the minimum 9 mm of the lower surface 
A_(0), i.e., 




FIG. 4: Diabatic potential surfaces V+ (dashed) and V- 
(solid) outside (left) and inside (right) the range where Zener 
flip tunneling occurs. 



that the tunneling rate follows, in case of vanishing dissi- 
pation, from the imaginary part of the partition function 



H = 



Pe 
2M 



V+{6) A(9) 



(53) 



Here, the diabatic potential surfaces read 



V± 



cos(6>) - i b 9 ± 




Vo[V(9)-V } 



1 o 



and the off diagonal elements 



A(9) = 



e c [Vo-V(0)} 



V' 2 



(54) 



(55) 



where V{9) = ejcos(^) and V = V{9 min ). By con- 
struction, for 9 — 9 m [ n the off-diagonal elements vanish. 
Further, one approximates these diabatic surfaces in the 
well-barrier range by a cubic expansion around the well 
minimum m ; n 

V±{q) = — ^ {q-q±) 2 [l-{q-q±)/q ,±}+6±,+ AV min , 

(56) 

where q = 9 — m i n is measured relatively to the mini- 
mum, 8^^ denotes the Kronecker symbol, and AV^n = 

V+(6*min) — V-(9 min ). 

Now, as long as V± are sufficiently separated from each 
other everywhere in the barrier range, i.e. \V— — V+\ 3> 
A, the spin degree of freedom is essentially frozen and 
the particle escapes through V- via standard MQT (see 
fig.0] left). The influence of noise for this scenario is thus 
captured by the findings of the previous sections. 

There is an additional domain, however . which gives 
rise to an interesting MQT phenomenon 19, 20] . Namely, 
when the two diabatic surfaces cross each other in the 
barrier range (see fig. El right), the particle's spin may 
flip during the tunneling since this may enhance the prob- 
ability for escape substantially. The standard MQT ap- 
proach has then to be extended to include spin flips so 



ft/3/2 



Z = Vqe 



s-[q] 



d,T\ 



-ft/3/2 



-hp/2 



x A(q(r 2 „)) • ■ • A(q(r 1 )) 

T2k 



x exp 



n '2k \ 

J2 J dr{V^q)-V + {q)) , (57) 



where S- denotes the bare action on the surface V- ■ The 
corresponding MQT rate can be cast into 



Ttot = T b + / flip (A c ) exp(-S flip /7i) (58) 



where Tt denotes the standard MQT rate without spin 
flip, while the second term captures contributions due 
to flips. The prefactor /flj p (A c ) with A c the coupling 
at the crossing point of V± is related to the probability 
for a spin flip to occur and Smp is the action along the 
flip-bounce. It turns out that a theory based on the non- 
dissipative partition function (|57fl provides rates which 
are already in good agreement with experimental data 
ppf. Nevertheless, a deeper understanding of the im- 
pact of low frequency noise on the Zcner-flip-tunncling is 
clearly needed. 

It was found in that when the particle traverses 
the intersection range of V± sufficiently fast the flip con- 
tribution is dominated by two spin-flips. Further, in case 
of weak friction the most profound effect of dissipation on 
the rate is provided by the action factor. Hence, we focus 
on the two- flip bounce and its action at low temperatures. 
In generalization of JSJ , the dissipative flip-bounce action 



takes the form 
Sfiip[g;si,S2] = 

hp/2 



'It iy? 2 + V-(q) + h Sl , S2 (r) [V + (q) - V-{q)\ 



-np/2 



hp/2 



dTd T 'q(T)K(T - r')g(r') (59) 



-hp/2 

where h Sl . S2 is the step function being unity in the inter- 
val [si, sa] and zero anywhere else and the kernel follows 
from (|39[) . For the two flip contribution the action also 
depends on the flip times si and S2- The influence kernel 
K(t) gives according to the analysis of Sec. IIV Al rise to 
a static and a dynamical contribution, see l|39fl • We first 
concentrate on weak friction for both, such that a per- 
turbation theory to first order in J s and k applies. By 
resorting again to a semiclassical evaluation this amounts 
to the fact to insert the undamped bounces gb and qfu p , 
respectively, into the influence functional. The rate en- 
hancement due to Zcncr flips can then be estimated by 
the difference between the simple bounce and the flip 
bounce actions, i.e., 



AS{J s ,k) 




= 0,K = 0) 



dTq b {r) 



Sfu p (J B 
hp/2 



= 0,K= 0) 



-i 2 



hp/2 



rfrgflip(r) 



(60) 



'-hP/2 



dTdr'k(T - t') [qh{r)qb(r') - qRi P (r)qm P (r')] 



The explicit form for the flip bounce and its action has 
been given in fl9l | to which we refer for further details. 
Since the flip bounce always exhibits are smaller ampli- 
tude and width (see fig-EJ), the terms in curly brackets are 
always positive. Hence, we arrive at the important result 
that weak static noise always lowers the rate enhance- 
ment due to Zener flip tunneling, while weak dynamical 
noise always increases its effect. The impact of which 
actually prevails depends only on the coupling strengths 
since the respective tunneling orbits enter both correc- 
tions in the same way. We note that so far (|6f ll relies 
only on the weak coupling limit, but applies to any kind 
of spectral density of the environment. 

As we have seen above, for a 1// noise the impact of 
dynamical modes is much smaller than that of the slug- 
gish ones so that we expect a shrinking of the enhance- 
ment. In order to see that explicitly, the tunneling orbits 
in presence of damping are determined numerically. Par- 
ticularly, the flip-bounce is determined by a vanishing 
variation of S[q; sj., Sa] with respect to all periodic paths 
and their corresponding flip times. It turns out that the 
optimal spin flips occur just as the trajectory reaches the 
intersection point and that they are centered symmetri- 
cally around r = 0. For the stationary action the time 




FIG. 5: 
(thick). 
(V+ thin 



Ordinary bounce (thin) and flip bounce trajectory 
The inset displays the inverted diabatic potentials 
, V- thick) intersecting at q c . 



=5 

< 




9/9 

FIG. 6: Rate enhancement due to Zener flip tunneling vs. 
temperature in presence of a bath with a 1/f spectrum and 
various values for the coupling constant (from top to bottom) 
j = 0.0 (solid), 0.1, 0.2, 0.3 (dashed). 



dependent potential can therefore be replaced by a piece- 
wise potential V c (q) = V^(q)Q(q c — q) + V + (q)Q(q — q c ) 
where q c denotes the intersection point. The equation of 
motion for the flip bounce thus reads 



- M ?ffip(r)+K'(«ffip) 



hp/2 



-hp/2 



dr' K{T-T')q mp (T') 



(61) 

Unfortunately the straightforward approach devel- 
oped for standard MQT rates and exploited also in 
Sec. IIIII where the equation of motion is solved numeri- 
cally by switching to Fourier space, does not work here 
due to the piecewise potential. Its derivative gives rise to 
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S function contributions, which pose numerical instabil- 
ities that cannot be circumvented by the usual rescal- 
ing. We thus determined qaip in coordinate space in 
an iterative procedure. One starts with the known flip- 
bounce for vanishing dissipation on the right hand 

ik) 

side and calculates the dissipative flip-bounce q^ by suc- 
cessively using <7gj p ^ for the dissipative part. The final 
flip-bounce follows after a sufficient number of iterations 
and by preserving the boundary conditions qfn p (T — > 
oo) = for HP — > oo. Typically, convergence is achieved 
quickly after 5 — 10 steps depending on the dissipation 
strength. From the flip-bounce the corresponding action 
is obtained which, for weak friction, leads together with 
the non-dissipative prefactor /fli P (A c ) and the standard 
MQT rate to the total rate in the flip range. 

In fig. IH1 the rate enhancement due to Zener flip tun- 
neling given by the action difference between standard 
bounce action and flip bounce action is shown for var- 
ious values of the coupling to a 1// bath. It is clearly 
seen that the impact of the static modes prevails and 
suppresses the enhancement with increasing temperature 
according to the above discussion. For weak coupling, 
however, the suppression remains small even at higher 
temperatures. To further illustrate the Zener flip effect 
on rate measurements, we show in fig.|7|the probability to 
escape P(ib) with varying bias current % [see ()30fl] for the 
quantronium circuit and in the range where Zener flips 
occur. Qualitatively, Zener flips lead to a pronounced 
shift of the s-curve towards smaller values of the exter- 
nal bias current, though, with a smaller slope in the first 
part of the curve. The influence of sluggish bath modes is 
to slightly move the Zcner-curve back towards the stan- 
dard MQT curve, but for weak damping this effect is 
hardly visible in the s-curves. By comparing with fig. [21 
however, one realizes that the Zener-curve resembles a 
standard MQT curve in presence of a stronger low fre- 
quency environment. Experimentally, some further in- 
formation about the influence of environmental modes 
is thus required to detect the Zener flip effect in a sin- 
gle s-curve measurement. The effect is clearly observable 
when one measures the bias current needed to maintain a 
certain value for P(ib) while sweeping through the Zener 
flip range by varying the magnetic flux <j> |2("|. 



VI. CONCLUSIONS 

Macroscopic quantum tunneling has regained new in- 
terest in the context of quantum information processing 
based on Josephson junction circuits. In this paper we 
analyzed the influence of low to moderate frequency noise 
on the tunneling rate out of a metastable well, where 
the noise spectrum is restricted to frequencies somewhat 
smaller than the typical frequency for oscillations around 
the well bottom (plasma frequency). A sluggish bath 
leads to a rate enhancement linear in the temperature, 
in contrast to the typical T 2 behavior for ohmic environ- 




FIG. 7: Probability to escape vs. external bias current for 
the quantronium circuit and in the parameter range where 
Zener flip tunneling occurs. The solid line is the result in- 
cluding Zener flips, while the dashed one is generated using 
the standard MQT rate. 



ments. For a sluggish bath with an additional 1// char- 
acteristic it turns out that even in presence of dynamical 
modes the static component by far prevails. This veri- 
fies that for MQT processes 1// noise can basically be 
treated as classical noise and in an adiabatic approxima- 
tion where for very weak coupling one effectively averages 
over rates for an ensemble of potential barriers. 

Based on these findings the impact of noise on a par- 
ticular kind of MQT process has been studied, where a 
particle carrying a spin-1/2 as an internal degree of free- 
dom experiences a Zener transition while tunneling. By 
comparing this new decay channel with ordinary MQT, 
we find that weak static noise suppresses the enhance- 
ment of Zener flip tunneling, while weak dynamical noise 
increases it. Further, we find that for escape probabili- 
ties (s-curves) the Zener flip effect leads qualitatively to 
a similar result as a standard MQT process under the 
stronger influence of sluggish bath modes. For the ex- 
perimental observation of Zener flip tunneling it is thus 
advantageous to have an environment weakly coupled to 
the system (higher Q-factor) and a spectrum with a de- 
creasing intensity towards low frequency modes, e.g. with 
a prevailing ohmic characteristic. 
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